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Q ■ Abstract 

h^ ! We show that gapped Graphene, with a local constraint that current arising 

^ I from the two valley ferniions are exactly equal, shows a non-BCS type super- 

. ' conductivity. Unlike the conventional mechanisms, this superconductivity 

c^ ■ phenomenon does not require any pairing. We estimate the critical tem- 

perature for superconducting-to-normal transition via Berezinskii-Kosterlitz- 
i-rt '. Thouless mechanism, and find that it is proportional to the gap. 
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1. Introduction 



(N 



Graphene is an atom thick allotrope of Carbon. Owing to its hexago- 



l/^ . nal honeycomb lattice structure it has unusual electronic properties. As was 

first shown by Wallace, that the valence band and conduction band touch 

^ I each other at six points in the momentum space, around which quasipar- 

ticle dispersion is linear [1]. It was pointed out by Semenoff that, the low 
energy quasiparticle excitations around these points, out of which only two 
are inequivalent, actually satisfy the massless Dirac equation [2]. Although, 

k> ', Lorentz symmetry is not a true space time symmetry to start with, it emerges 

H I in the long wavelength limit because of honeycomb structure of the material, 

the velocity of light being replaced by Fermi velocity of the electrons. Ex- 
perimental isolation of this material paved the way to simulate relativistic 
quantum physics and check validity of various fertile ideas [3]. The phe- 
nomena that make this material important include Klein paradox [4], room 
temperature quantum Hall effect [5], Andreev reflection [6] and universal 
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conductance [7, 8]. Fermion number fractionalization due to vortices [9] and 
anyonic excitations [10] have also been studied in this novel material. 

Although, the experimental search for superconductivity in pure graphene 
is yet to be successful, possibility of superconductivity in graphene has at- 
tracted considerable attention from both theoretical and experimental sides 
[11, 12]. Experimental observation of proximity induced superconductivity 
in graphene, implies that graphene supports phase coherence [13]. Variants 
of conventional pairing mechanisms are proposed for superconductivity in 
pure and doped graphene, and graphitic layers [14, 15]. In Ref. [16], the 
authors address the possibility of high Tc superconductivity using resonating 
valence bond model for doped graphene. Doped graphite is known to be su- 
perconducting, and possible mechanisms giving rise to this effect have been 
proposed [17, 18]. In Ref. [19], the authors propose a mechanism of super- 
conductivity arising due the edge states in graphene. In Ref. [20] authors 
address a possibility of non-pairing based superconductivity. 

In the present paper, we show that Graphene with a finite non-zero band 
gap and a constraint that, the currents from both the valley fermions are 
always equal, exhibits superconductivity. We find that Graphene with this 
constraint, possesses infinite DC conductivity, shows Meissner effect and flux 
quantisation. It is seen that, the Lagrange multiplier field introduced to im- 
plement the above constraint behaves like a Nambu-Goldstone mode of BCS 
theory, and plays the central role in realising superconductivity. However, 
unlike BCS theory, here Meissner effect and flux quantisation occur not due 
to Anderson-Higgs mechanism, but rather due to topological Chern-Simons 
coupling. We find that the full quantum theory has no propagating Dirac 
fermion, and only charge neutral fermion-hole bound pairs propagate. After 
a certain finite temperature, we observe that spontaneous proliferation of 
monopoles in Lagrange multiplier field takes place via Berezinskii-Kosterlitz- 
Thouless phase transition, which marks the superconductor-to-normal tran- 
sition. In case of finite size Graphene with armchair egdes, we show existence 
of dissipationless chiral gapless edge modes, as a consequence of supercon- 
ductivity in bulk, . 

The paper is organised as follows. In subsequent section, we discuss in de- 
tail both classical and quantum theory of Graphene with current constraint. 
A study of electromagnetic response of this theory is done in section 3. Sec- 
tion 4 deals with discussions about phase transition in this theory, followed 
by section 5 where boundary theory is worked out. Section 6 summarises the 
results and closes with a discussion about obtained results. 



2. Current constrained model 

In tight binding approximation, the electron hopping Hamiltonian, de- 
fined on a hexagonal lattice with hopping energy t, reads: 

H=-J2Y1 MaHr)6(r + Si) + 6t(r + Si)a(r)), (1) 

r i=l,2,3 

where the fermion operators a and b act on the two interpenetrating sublat- 
tices A and B ^. As is well known, H can be linearized around two Dirac 
points K± ^, to yield: 



Hd = -ivF / d'^x |^+(r)o- ■ V^+(r) 



+^Ur)ff*- V^_(rl 



where vp = ^ is the Fermi velocity (which is hence forth set to unity, 
along with h), ^E'±(r) = (a!t(r), 6j_(r)) and Pauli matrices are defined as 
a = {(Xx, ay) Sz a* = (cr^, —cfy) [1, 2]. By selective doping of one sublattice or 
due to interaction with substrate, a local onsite potential can be generated 
which breaks sublattice symmetry, and leads to a mass gap in the electronic 
spectrum. This renders Graphene a semiconducting behaviour, and if mass 
gap is large enough then insulating nature follows. The Lagrangian describing 
these massive Dirac fermions can be written in a manifest Lorentz invariant 
form as: 

^D = ^+{h+d^ - m)ij+ + V5_(?7^(9^ - m)ij_. (2) 

Here, we have defined Gamma matrices for ■?/'+ field as 7° = o'3,7+ = 
iai and 7^ = ia2- Gamma matrices for ■?/'_ field are also same as ■?/'+ ex- 
cept for 7^, which is defined as 7^ = —7^ (henceforth we shall use Feynman 
slash notation, whereby (^^ = 7±a^). So we see that in Graphene, low energy 
electronic excitations are two species of Dirac fermions concentrated around 
two {K±) valleys. We now consider a scenario where these quasiparticles are 
coupled to an interaction, which leads to preservation of valley symmetry 



^Hcre vector r points to sublattice A, which is connected to B sublattice via vectors 



si = (1, V3)|, S2 = (1, -\/3)^ and S3 = {-1,0)1, with I being the C-C bond length. 

2tt _|_ 2tt ^ 
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Hn our convention K± = (|2:,±-2^). 



locally. In what follows, we shall not be concerned about the detailed struc- 
ture of this interaction, but make use of fact that, that currents generated 
in response to some external stimuli, from both valleys must exactly be the 
same locally i.e., j^{x,t) = jt{x,t) {i^{x) = ■ilj{x)'j^ip{^))- Further, we shall 
also assume that, although the fermion fields are strictly interacting, the in- 
teraction allows above mentioned Dirac description of quasiparticles, and the 
net effect of interaction amounts only to the above local current constraint 
(LCC). 

The above condition endows the theory with a local gauge invariance. 
Indeed it is straightforward to check that, under following local gauge trans- 
formations: 

^+ ^ e-^^(^V+ (3) 

ip. -^ ilj_e'^^''\ (4) 

the Lagrangian transforms as: 

^^^^^ + (j^l-j^)d,e, (5) 

and since j1 — jt = 0, the Lagrangian remains invariant under local gauge 
transformation. Hence, we have shown (at the classical level) that, the theory 
under consideration is an Abelian gauge theory, albeit without gauge field 
[21]. 

In functional integral formulation of quantum field theory, generating 
functional is an object of central importance, which for the above mentioned 
theory reads: 

Z = nJ ^[V5(+,_), ^(+,„)] 6 {f+ - f-) e''-, (6) 

where Sd = I d^x^n- (7) 



Delta function in the above expression is introduced to implement LCC, and 
can be rewritten by introducing an additional Lagrange multiplier Bose field 



a^, such that: 



Z[^(+,_),^(+,_)] = N / ^[V5(+,_),V^(+,_),a^]e^^-[^'^''^-], (8) 



where 

+^lj_{i$_-m-(fiJi)J\. (9) 

Modulo an unimportant normalisation constant, (6) and (8) describe same 
physics, as can be checked by integrating out a^ field to reproduce the delta 
function constraint. This shows that above action captures dynamics of 
quasiparticles in Graphene subject to LCC. It is evident from the structure of 
the above action, that it remains invariant under local gauge transformations: 

^^ ^ e-^^(^V+, i'- -> ^^e^'^"\ a^^a^ + d^d{x). (10) 

Hence, one can identify a^ as an Abelian gauge field, which is minimally 
coupled to two fermions, which are oppositely charged. It is easy to see 
that generating functional defined above is also invariant under above trans- 
formations, and no non-trivial Jacobian appears since these transformations 
are assumed to be regular. As a consequence, we expect that various n-point 
functions in this theory would obey Ward-Takahashi identities. In our case, 
two-point function for i^+ valley fermions, is given by: 

iSpix -y) = {T {i)+{x)i}+{y))) 

= N [^[ij+,ij+,a,] ^+(a;)^+(i/)e*^-['^''^'«-l. (11) 



Under field redefinition (10), we see that propagator satisfies a Ward-Takahashi 
identity, 

Sp{x -y) = e-''^^^^SF{x - y)e''^y\ (12) 

whose only solution is Sp^x — y) oc 5{x — y). Above identity is very powerful, 
since it has allowed for an exact determination of propagator in this inter- 
acting theory. Exactly similar identity would also hold for propagator of K_ 
fermions. It is worth mentioning, that this model is one of the rare cases 
where full propagator of this theory is known without any approximation. 
Presence of a physically observable particle in a theory, manifests as poles 
of propagator in momentum space. In our case, as is clearly evident, the 
propagator is regular everywhere in momentum space, which implies that 
Dirac fermion in our theory is not a propagating mode. This is particularly 



surprising since we started with a free Dirac theory with a constraint condi- 
tion on currents, and it appears that condition is severe enough to not allow 
free fermion propagation. 

In the absence of Dirac fermions, it is a natural to inquire about quasi- 
particle excitations in this theory. Inorder to answer this question, it is 
instructive to study the four-point function in this theory, which is defined 
as: 

(T {ilj+{xi)ilj+{x2)ip+{yi)ip+{y2))) = 

N j ^[^+,^+,a,] V^+(xi)V^+(a;2)V5+(i/i)^+(|/2)e^''". (13) 

Under local gauge transformation (10), we obtain another Ward-Takahashi 
identity for four-point function: 

{T {i)+{xi)i)+{x2)i)+{yi)i)+{y2))) = 

e-^'(^^) g-.^{..)^2- (V;+(a;i)V^+(x2)V+(l/i)V^+(y2))) e^'^^^^ e''^y^\ (14) 

Apart from a trivial non-propagating solution discussed above, assuming 
validity of translational invariance, above equality admits a solution of the 
kind: (T (^+(a;i)^+(x2)^+(i/i)^+(|/2))) oc 5{xi - yi) 5{x2 - 1/2) /(a^i - X2), 
where / is some function of {xi — X2). This means that the above identity 
allows for propagation of composite operator il){x)il){y)\x=y, which describes 
charge neutral excitations consisting of fermion-antifermion bound states. 
In this case, these can be conveniently identified with exciton excitations 
of Graphene. Hence, we see that instead of Dirac fermion, Graphene with 
above constraint, admits excitons as its quasiparticle excitation. It is worth 
mentioning, that absence of fermions as elementary excitations and occurence 
of bound states in a constrainted theory like above, also appeared in a model 
of color confinement proposed by Rajasekaran and Srinivasan, and in related 
works [21, 22, 23]. Interestingly, they showed that quarks and gluons (which 
appeared as bound states) did not propagate and were confined, whereas 
mesons (color neutral bound states of quarks) were propagating excitation 
in their model. 

3. Electromagnetic response 

The external photon field, living in 3 + 1 D space time, interacting with 
fermions confined on a plane can not be correctly described by —^F^^F^'^ 
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(/i, z/ = 0, 1, 2), since it has only one degree of freedom, and not all compo- 
nents of physical photon field couple to matter. As shown in Ref. [24], one 
begins with action for electromagnetic field j d^x^F^^F^^'^ + j^A'^, where 
j^ describes the matter current confined to XY-plane. Integrating out z- 
coordinate in the above action, and using Greens function identity, we get 
the three dimensional action as : / dP'x^F^y^=F'^'' + j^A^, where F^^, (/i, 
z/ = 0, 1, 2) describes 3D physical magnetic field Bz-, in-plane components of 
electric field E^ and By. So, Lagrangian (9) with suitable coupling to photon 
field reads: 

^ = - A^M.^^^'^ + [^+(^^+ -m + ^ + 4)^+ 



+ilj_{i$_-m-^ + 4)'ip-]- (15) 

It needs to be pointed out here that, we have assumed velocity of light c equal 
to the Fermi velocity vp which is set to unity. This is actually not the case, 
but as shown in Ref. [25, 24], this does not change the qualitative behaviour 
of the theory. Also, notice that the electromagnetic coupling constant g is 
not dimensionless, as in QBD4, but rather has dimension of square root 
of mass. Above Lagrangian describes our theory at tree (classical) level. 
Inorder to take into account effects due to quantum corrections, which arise 
from virtual fermion loop excitation, one needs to find out the effective action 
by integrating out fermion field. Fermion spectrum is gapped in our theory, 
and at low energies fermions are only excited virtually, hence integration of 
fermion fields is physically meaningful. Effective action upto quadratic terms 
in fields, obtained using derivative expansion of fermion determinant [26, 27] 
reads: 

67r|m| T[\m\ 

where an additional factor of 2 has been multiplied to account for spin de- 
generacy of fermions. It can be shown that, in the limit of large m this 
approximation is valid and higher order terms can be neglected. As is evi- 
dent, a^ did not have a kinetic term to start with, but fermion loops have 
made it dynamical, and indeed it can be identified as a genuine Abelian gauge 
field. Further, these two gauge fields are coupled by a mixed Chern-Simons 
term, which has a topological nature [28, 29, 30, 31, 32]. In other words, 
this implies that a^ field has now become electromagnetically charged due 
to presence of virtual fermion cloud around it, with current being given by 



J^ = t^'^f'dyap. It is interesting to note that, by defining A± = A ± a fields, 
the Chern-Simons term in above Lagrangian can be rewritten, modulo a con- 
stant, as : e^upA^d" A^j^ — e^upA'^d^A'L, which has two topological terms with 
opposite signs; as a result it does not violate parity. Hence, the net action 
describing interaction of Graphene with electromagnetic field is given by: 

S=/d»x^, where 



^ = -\f^-^F"" - T-^-iM"" - -T^e"-M,a„a,. (17) 

4 ^(}2 67r|m| 7r|m| 

Note that, above action is invariant under two types of local gauge trans- 
formation: a^^ ^ a^ + d^\ and A^ — )> A^ + d^Xi where A and x ^-^e some 
regular functions of x, and we have assumed that fields under discussion de- 
cay sufficiently quickly as one approaches Graphene boundary. In general, 
it is easy to show, that these gauge invariance would imply presence of a 
massless mode living on Graphene boundary. However, in what follows, we 
shall assume that Graphene sheet is practically infinite, and discussion about 
bulk-boundary interaction will be pursued later. Further, first two terms in 
above Lagrangian are invariant under these transformations, whereas the 
third Chern-Simons term is not. Inorder to appreciate consequences aris- 
ing due to this noninvariance, we follow Ref. [25] and do a Wick rotation: 
t — 7- —IT (where r G [0,/9] is a compact variable), from Minkowskii space- 
time to Euclidean space-time. In this compact Euclidean space-time, bosonic 
(fermionic) fields obey periodic (anti-periodic) boundary conditions: 

Bosonic: B{x, 0) = B{x, P), (18) 

Fermionic: ■?/'(x, 0) = —■?/'(x,/3). (19) 

These conditions alongwith requirement of single valuedness, impose a re- 
striction on gauge function \{x): X{/3) = A(0) + 27in, (n being integer, often 
called the winding number). Vacuum functional is defined as: 

^Euclid _^ f ^A^^a^ e-^c5 , where (20) 

Scs = f dr f £x ^e^-PA^d,ap. (21) 

Jo J 7r|m| 



Considering variation of Chern-Simons action under restricted gauge trans- 
formations where A(r) only depends on r, one finds: 

SScs = -T^n^, where (22) 

\m\ 

^= f £xe,jdiAj, (23) 

is magnetic flux. Invariance of vacuum functional (20) under above transfor- 
mation, demands: SScs = '^'^'iN, where N is an integer. This clearly implies 
that $ = ^nUil which when written in SI units reads: 

$ = iV ^ ^ — . 24) 

\\m\)2g ^ > 

Hence, magnetic flux is quantised in this model with flux unit |2. 

Inorder to understand response of the system due to influence of external 
electromagnetic fleld, we integrate out a^ fleld from Lagrangian (15), to arrive 
at an effective action for electromagnetic fleld: 

3|m| /, ,„ , (9^9^^ , \ 
^eff = -^ [A,A^ - A.^A^j . (25) 

As is evident, interaction with a^ fleld has induced mass M = -^ for the 
physical electromagnetic fleld. The photon fleld becoming massive implies 
Meissner effect, a hallmark of superconductivity, where the static magnetic 
fleld exponentially dies down with distance from the boundary, with charac- 
teristic length scale, called penetration depth, which in our case is A = g^. 
Further, current- current correlation function in this theory, can be found 
from above Lagrangian to be: 



{f{x)f{y)) 



5^S, 



eff 



5A^{x)5A,{:y) 



3|m| / „„ d^^d^ \ ^n, , , ^ 

4^(^" -^^)^'(^-^)- (26) 



Notice, that above function has a pole at zero momentum, which implies 
inflnite DC conductivity from Kubo formula. To make the comparison with 
usual BCS theory clear, we follow the technique given in Ref. [33], and 
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define an auxiliary scalar field 0, such that d^cj) = e^^pd'^af, so that the 
above Lagrangian after a Hubbard-Stratonovich transformation, reads: 

37r|m| /„ , m _, \ ,^„, 

Inorder to make comparison with Landau-Ginzburg theory apparent, one can 
define x = 6**^, so that the above Lagrangian reads: 



_ SttH 

■^eff - 7^ 



2iT\m\ 



(28) 



As is clear, the above Lagrangian is in manifest London form, and all 
the phenomenological properties of superconductivity would follow from here 
[34]. It was shown in generality by Weinberg [34], that occurrence of an elec- 
tromagnetically charged field ^(x), which transforms as ^{x) — )■ ^(a;) + A(a;) 
under a gauge transformation, is sufficient for the theory to exhibit super- 
conductivity. In case of theories which exhibit any kind of fermion pairing, 
phenomenon of spontaneous symmetry breaking takes place, whereby system 
realises a vacuum, which is not invariant under all symmetry group transfor- 
mations that preserve the Hamiltonian. If the symmetry not respected by 
vacuum, is local gauge invariance of electrodynamics, then there inevitably 
appears a massless Nambu-Goldstone field, which transforms exactly like 
^{x) — )■ ^{x) + A{x) under a gauge transformation [35, 36]. As argued by 
Weinberg, it is presence of this Goldstone mode that is ultimately responsi- 
ble for superconductivity. One can see from above Lagrangian, that in our 
theory, field behaves exactly like the Nambu-Goldstone mode, and hence 
our theory exhibits genuine superconductivity. However, it is worth noting, 
that in our theory, the origin of field is due to a local constraint, and not 
due to breaking of any local symmetry by vacuum. In other words, our theory 
exhibits superconductivity due to presence of a local constraint, rather than 
that of fermion pairing. Also it is worth pointing out that, this non-BCS 
type superconductivity relies on presence of Chern-Simons term and hence 
is present only in the planar 2 + 1 dimensional world. 

4. Phase transition 

In above discussion, we have tactically assumed that a^ field is regular and 
single valued everywhere. In general, this may not be the case, which allows 
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for occurrence of monopoles corresponding to a^ field, i.e., e^ypd^f^^ ^ 0. 
Equivalently this translates to presence of vortices in field, which means: 

<h V0 ■ ds = ±2nn, (n being an integer) (29) 

along any closed curve C that enclosed the vortex. Interestingly, in absence 
of Afj^ field, above action exactly matches with that of continuum limit of 2D 
XY model. In case of XY model, it was shown by Berezinskii [37] and inde- 
pendently by Kosterlitz and Thouless [38] that, presence of vortex excitation 
is forbidden energetically at zero temperature, since it does not correspond 
to minimum of free energy. At zero temperature, vortices and antivortices 
are bound together due to mutual attractive potential, and hence do not 
move freely. However, at temperatures above a certain critical temperature, 
they showed that, presence of vortices is energetically favourable as it min- 
imizes free energy. This means that at critical temperature, a transition 
from bound vortex state to free vortex state takes place, which is known as 
Berezinskii-Kosterlitz-Thouless (BKT) phase transition. This is an infinite 
order phase transition, wherein the ordered(low temperature) phase differs 
from disordered(high temperature) phase in terms of topology of field, rather 
than symmetry. It was shown that presence of these vortices, destroys the 
long range correlation present otherwise in the system. It is also known that 
the ordered phase in case of XY model, does not exhibit any Long Range 
Order, but shows what is called a Quasi Long Range Order, whereby the 
correlation function exhibits a power law decay asymptotically. So as in case 
of XY model, in this case also one expects a BKT phase transition to take 
place, after which free vortex excitations could be present. The critical tem- 
perature for this transition, can be read off from above action, which is given 
by 

Tbkt = -^-^. (30) 

Note, that the critical temperature depends linearly on mass gap, and for 
gapless fermionic spectrum, the transition would occur at zero temperature. 
The effect of vortices on electromagnetic response can be understood by 
writing an effective action for electromagnetic field after taking into con- 
sideration occurrence of vortices. Following Ref. [39, 40], we take vortex 
excitation into account, in above action by writing — )■ (p^eg + 4>vor, where 
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(preg{vor) IS regular (vortcx) part of 0, so that action reads: 

■^eff = ^ ( df,(j)reg + df,(l)yor + Y~\ 1 ^M J ' (^l) 

Using an auxiliary field C,^, above can be rewritten as: 

Integrating out regular part of 0, implies a constraint: d'^C.fi = 0, which has 
an obvious solution ^^ = e^^pd'^B^, where B^ is some gauge field. Above 
action written in terms of B field is given by: 

+ B, U^'d,d,(j^,or + '^^^^e^^^d,A)j . (33) 

Integrating out auxiliary B field, one gets an effective action describing in- 
teraction of vortex current K^^ = e^^''dudp(f)vor with electromagnetic field A^: 



37r|m| 



4 









(34) 



Notice, that the last term in above expression, is topologically trivial and 
is responsible for superconductivity. In absence of vortex current, it solely 
describes the response of this system. First term on the otherhand, describes 
interaction between vortices, and it can be easily shown that interaction po- 
tential between them has a logarithmic behaviour. The other two terms, in 
the effective action, describe interaction of vortices with electromagnetic field, 
which means that vortices in our model are charged. So, in order to find con- 
tribution of these charged vortices to electromagnetic response, one requires 
to integrate out vortex field, to yield a net effective action for electromagnetic 
field. Very interestingly, integration over vortex current contributes a term 
of the form A^ (^^ — i]^'^) A^ which exactly cancels with the one already 
present in above action. This straightforwardly means that, in presence of 
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vortices, superconductivity is destroyed [40]. Hence, in this model, we see 
that occurrence of BKT phase transition marks a superconducting-to-normal 
phase transition, and therefore Graphene with local current constraint, at 
sufficiently low temperatures, exhibits superconductivity with strong type-II 
character. 

5. Boundary theory 

As mentioned earlier, we have assumed that Graphene sheet is of infinite 
extent, and fields fall of sufficiently quickly, so that surface terms give negligi- 
ble contribution. However, that is rarely the case in reality, where in general 
one encounters Graphene samples with boundary. Owing to its hexagonal 
tiling, Graphene can exhibit boundary of two kinds: Arm chair and Zig-zag. 
It is known that, the latter exhibits localised electronic egde states, whereas 
the former does not. Hence, in case of arm-chair egdes, fermions present in 
bulk can freely interact with the ones living on boundary and vice versa. In 
what follows, we shall assume that Graphene sheet has a well defined regular 
boundary of arm-chair kind. 

As noted above, low energy effective action describing dynamics of low 
energy electronic excitation, subject to the local constraint, coupled to elec- 
tromagnetic field is given by: 

1 TTl 
"^ = ~7 — 1 iJuuJ j fC ^nOudp- 

b'n\m\ Tx\m\ 

As was observed earlier, the last term in above Lagrangian is not invariant 
under local gauge transformation: a^ — )> a^ -|- 9^A, where A is some regular 
function of x. As a result, the change in action is given by: 

5Scs = (^^) / d'x e^^^d, (A /.,) . (35) 

Above volume integral can be converted to a surface integral, defined on 
closed Graphene boundary, to give an action: 



55, 



cs 



(^)//xe-A/„„, (36) 



This term, as it stands, is not gauge invariant, and is defined on Graphene 
boundary, which encloses the bulk. Gauge invariance of any given theory, is a 
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statement that, the theory is constrained, and possesses redundant variables. 
We observe that, our theory to start with was gauge invariant at classical 
level. One loop corrections arising out of fermion loops, generate Chern- 
Simons term, which exhibits gauge noninvariance. Because, our theory to 
start with was gauge invariant, and hence constrained, consistency demands 
that quantum(corrected) theory should also respect the imposed constraints, 
and hence should be gauge invariant. The occurrence of above gauge nonin- 
variance, simply implies that one is only looking at one particular sector of 
theory, and there exists other dynamical sector, whose dynamics is such that 
it compensates with the one above to render the total theory gauge invari- 
ance. Following Ref. [41], we demand that there must exist a corresponding 
gauge theory living on the boundary, defined such that it contributes a gauge 
noninvariant term of exactly opposite character and hence cancels the one 
written above. The simplest term, living on boundary, that obeys above 
condition is: 

Sb = ^^^^ [ d'xOe'^'^f,., (37) 

^TT Jb 

where 6{x, t) is Stiickelberg field, which transforms like 6 — ?> 6+ A. In general, 
this scalar field would be dynamical, and with a gauge invariant kinetic term, 
the boundary action reads: 



'B 



Sf! = I (fx 



Zn 



(38) 



Note, that because of its peculiar transformation property, a quadratic mass 
term for 6 is not gauge invariant. Hence, in a gauge theory framework like 
this, 6 field remains massless. In deriving above action, we have only consid- 
ered gauge invariance with respect to transformation in a^ field. However, 
analogously the same may be done for A^ field, so that net action, describing 
massless surface modes, coupled to both gauge fields is given by: 

Sb= f d^x [c {d^9 -a^- A^f (39) 

Jb 



sgn{m 



,^—ee^^ (/,. + F,^) 



(40) 



There are several things to note here. Firstly, the action for 9 is in mani- 
fest London form, and hence is indicative of non-dissipative transport on the 
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Figure 1: Figure depicting motion of K^i_\ valley fermions in red (blue), in Graphene 
under LCC. LCC is maintained at any point P in bulk, as it gets equal contribution 
from both the valleys in all directions. Points on the boundary, however, like point Q, 
get contribution only from Graphene interior, which leads to existence of clockwise and 
counter clockwise edge currents (shown by red and blue arrows), one from each valley, such 
that the net edge current adds up to zero. But, Graphene edge being one dimensional, 
possesses chiral anomaly, which gives rise to a net non-zero chiral edge current (shown by 
black arrow). 



boundary. Secondly, 6 field is electromagnetically charged, and hence bound- 
ary supports dissipationless electric current, or in other words boundary is 
superconducting. Thirdly, the coupling of 6 field, with that of a and A field 
is anomalous, as a result chiral current in this quantum theory is no longer 
conserved. This ultimately results in chirality of these surface modes. This 
can be explicitly seen by observing that, unlike bulk fermionic coupling, these 
surface modes couple to sum of two gauge fields a + A, and not to difference 
a — A. Hence, we have shown that, arm chair edged Graphene with LCC 
possesses dissipationless chiral gapless surface modes. 
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6. Conclusion 

We have shown that, gapped Graphene with a local constraint j^ = jt, 
exhibits superconductivity, albeit without any pairing. Apart from infinite 
DC conductivity, we show that this theory also realises Meissner effect, a hall- 
mark of superconductivity. However it is worth noting that, unlike BCS the- 
ory, here Meissner effect does not occur due to Anderson-Higgs mechanism, 
whereby photon becomes massive by absorbing gap less Nambu-Goldstone 
mode. Rather, the origin of Meissner effect (or presence of photon mass) 
is due to presence of topological Chern-Simons term in the theory. As 
is known, Chern-Simons term is unique to two spatial dimensions, which 
means that this type of superconductivity only occurs in planar world and 
its higher dimensional extension may not be possible. It is seen that, the 
Lagrange multiplier field introduced to implement above local constraint, is 
found to behave analogous to Nambu-Goldstone mode of BCS theory. In 
pairing based theories, because of spontaneous breaking of symmetry there 
appears an excitation mode called the Amplitude mode [42]. In our theory, 
there is no fermion pairing, which implies absence of amplitude mode. This 
feature makes this superconductivity distinguishable from the other pairing 
based theories. At a critical nonzero temperature, like in 2D XY model, we 
see that spontaneous proliferation of monopoles in Lagrange multiplier field 
takes place, and marks occurrence of BKT phase transition. It is seen that 
the contribution due to these monopole (singular) part of this field, exactly 
cancels the same coming from the regular part. So loss of superconductivity 
takes place via BKT mechanism, and superconductor-to-normal transition 
takes place at BKT critical temperature, which is proportional to mass gap, 
and hence can be controlled. On a finite graphene sheet with armchair egdes, 
we show that, this type of superconductivity naturally supports gapless chiral 
edge modes. Amusingly, we find that the full quantum theory does not have 
Dirac fermions as propagating mode, and charge neutral fermion-hole bound 
pairs (excitons) are found to be the elementary quasiparticle excitation of 
the theory. 

Above discussed counterintutive features of our theory can be under- 
stood by noting that the Lagrange multiplier field introduced to implement 
current constraint (equation (8)), minimally couples to Dirac fermions and 
hence manifests like an effective magnetic field. However, the sign of this 
magnetic field is opposite for both the valley fermions, which now move in 
mutually opposite cyclotron orbits in accordance to effective magnetic field 
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felt by them. As shown in Fig. 1, the current constraint is obeyed at any 
point in bulk, as the contribution to local current coming from cyclotron 
loops from both valley fermions would be equal. At Graphene boundary, the 
contribution from cyclotron loops present in vicinity adds up, giving rise to 
clockwise and counter clockwise circulating edge currents of equal magnitude. 
Classically(at tree level), these two currents add up to zero and hence no edge 
current exist. However, in one dimensional theories like Graphene boundary, 
quantum corrections generate terms at one loop level, which violate conser- 
vation of chiral currents. As a result of this, the two opposite current flows 
no longer cancel each other and a net nonzero edge contribution survives. 
Fermions at the two valleys living in bulk feel exactly opposite effective mag- 
netic field, so presence of an external magnetic field creates an asymmetry 
between the two valley fermions, which now feel different magnetic fields. 
It is this difference of magnetic fields that results in difference between cy- 
clotron orbits amongst two valleys, and leads to non zero circulating current 
in the bulk. This ultimately gives rise to net diamagnetic response of the 
system. A neutral composite particle like exciton would not couple to any 
of these fields, and hence would be a long lived quasiparticle excitation in 
absence of Dirac fermion, which does not show up in the Hilbert space of 
asymptotic fields. 

Above theory is reminiscent of a model of superconductivity, proposed by 
Laughlin [43], which goes by the name of anyon superconductivity [44, 45, 46, 
47] (for a lucid introduction and broad review see [41]). It was shown in this 
context that, fermions interacting with an emergent statistical gauge field 
can become anyons, and show superconductivity, albeit without any pairing. 
A crucial feature in which our theory differs from anyon superconductivity 
is that, in the latter case time reversal symmetry is explicitly broken due 
to presence of Chern-Simons term in action, which is responsible for frac- 
tional (anyonic) spin and statistics of quasiparticles. On the otherhand, in 
the present proposal, time reversal symmetry is preserved and quasiparticles 
do not have anomalous spin/statistics. Although above theory has strong 
resemblance with the one proposed in Ref. [20], it needs to be pointed out, 
that these two theories are fundamentally different, in the sense that in the 
latter case one does have Dirac fermion as a propagating mode, whereas in 
present case one does not. 
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